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Abstract 

This paper introduces a generalization of the so-called space-fractional Poisson process by extending the dif¬ 
ference operator acting on state space present in the associated difference-differential equations to a much 
more general form. It turns out that this generalization can be put in relation to a specific subordination of a 
homogeneous Poisson process by means of a subordinator for which it is possible to express the characteriz¬ 
ing Levy measure explicitly. Moreover, the law of this subordinator solves a one-sided first-order differential 
equation in which a particular convolution-type integral operator appears, called Prabhakar derivative. In 
the last section of the paper, a similar model is introduced in which the Prabhakar derivative also acts in 
time. In this case, too, the probability generating function of the corresponding process and the probability 
distribution are determined. 
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1 Introduction and background 


In the last decade, it became apparent that several phenomena that can be modeled in terms of point pro¬ 
cesses are non Poissonian in nature (see as examples). As a consequence, there has been an increased 

interest in generalizing the Poisson process iV(t). The Poisson process is a counting process with many nice 
properties. It is a Levy process and, therefore, its increments are time-homogeneous and independent. It is 
a renewal process, meaning that the sojourn times between points (or events) are independent and identi¬ 
cally distributed following the exponential distribution. It is a birth-death Markov process and its counting 
probability obeys the forward Kolmogorov equation 

d 

-:rPkM =-^PkM + ^Pk-iM, ( 1 - 1 ) 

dt 

where = P{iV(t) = fc}, fc > 0, t > 0 are the state probabilities of the Poisson process and A is the rate 
of the Poisson process. There are many possible ways to generalize this process. We are interested in the 
so-called fractional generalizations of the Poisson process where either the derivative on the left-hand side or 
the difference equation on the right-hand side of (11.11) are replaced by suitable fractional operators. 

For instance, if one keeps the renewal property and considers sojourn times such that ¥{Np[t) = 0} = 
for 0 < /3 < 1, where Ep[z) = /E(.kp + 1) is the one-parameter Mittag-Leffler function, 

one gets the renewal fractional Poisson process Np{t) discussed in This leads to the equation for 

p^it)=F{Npit) = k} 0 ], 

d^ 

■^PkW =-^PkM + ^Pk-iM, k>0,t>0, ( 1 . 2 ) 

(with A = 1) where d^/dt^ is the so-called Caputo derivative, a pseudo-differential operator defined as 


dl^ 1 

r(i-/3) Jo 


(t-u)- 


i d/(u) 

du 


du. 


The renewal fractional Poisson process is not a Levy process 


(1.3) 
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Another possibility is generalizing the Poisson process via the so-called space-fractional Poisson process 
studied in ll33ll . The space-fractional Poisson process is in practice a homogeneous Poisson process subordi¬ 
nated to an independent stable subordinator. We describe this model in details in Section ITT] Here we aim at 
a further generalization of this process preserving the Levy property. This generalized process is constructed 
by means of a superposition of suitably weighted independent space-fractional Poisson processes. The result¬ 
ing process is then subordinated by means of a random time process in order to account for the modelization 
of a possible irregular flow of time. The complete construction of the model is done in Section ITT 

The motivation at the basis of this study lies in the importance that superposition of, possibly dependent, 
point processes plays in applied sciences. An important example concerns the modelling of neurons’ incoming 
signals. It is commonly accepted that each neuron obtains information from the neighbouring neurons in the 
form of spike trains (i.e. signals with powerful bursts), closely resembling realizations of stochastic point 
processes. Therefore, each neuron, receives a superposition of, possibly rescaled, point processes and its 
subsequent behaviour depends on the characteristics of this afferent combined input signal. Based on classical 
result contained in H, 3] (roughly saying that a superposition of sufficiently sparse independent 

point processes converges to a Poisson process), the inpu t sig nal was considered in this applied literature to 
be well approximated by a Poisson process [see e.g. 1201.13811 . However, it was later shown that the above 
result does not always apply to superposition of signals from neurons’ activity and that experimental evidence 
deviates from a Poissonian structure [|2^. [lih . Recently, the study of weak convergence of superposition 
of point processes has regained interest [see e.g. |T3, and the references therein] showing that different 
behaviours are possible. Within this framework, we can consider the model we are going to describe as a 
weighted finite superposition of independent space-fractional Poisson processes (each of them generalizing 
the homogeneous Poisson process but also admitting the possibility of jumps of any integer order). Each of 
these space-fractional Poisson processes can be thought to model different groups of neurons (different areas 
of the brain) acting together with simultaneous spikes giving rise to the non-unitary jumps. Notice, finally, 
that the space-fractional Poisson process is a non-renewal process and so is the generalized space-fractional 
Poisson process. This is a key feature for the combined neurons’ input signal as it is remarked in ll2^ . 

For the sake of clarity and simplicity, we first recall some selected basic mathematical results regarding sub- 
ordinators which will be useful in the following. Section[2]presents the construction of a random time-change 
by means of independent subordinators and tempered subordinators. This is of fundamental importance for 
the definition of the generalized space-fractional Poisson process which will be carried out in Section[3] 

We start thus by recal ling some basic facts on subordinators and tempered stable subordinators. The 
reader can refer to or 02611 for a more in-depth explanation. We recall that a subordinator is an increasing 
Levy process defined as follows. 

Definition 1.1. Given a filtered probability space where ^^e associated right- 

continuous filtration, the process S^, t > 0, adapted to jj, starting from zero, and with increasing paths, is 
called a subordinator if it has independent and stationary increments or, equivalently. Vs > t > 0, — Sj is 

independent of and 8^—5^= in distribution. 

In this paper we consider only strict subordinators that is those with infinite lifetime and with the only 
cemetery state (oo) located at t = oo. For this special class of subordinators a simplified version of the well- 
known Levy-Khintchine formula for general Levy processes holds. In particular, the following theorem gives 
a characterization of subordinators in terms of Laplace transforms of their one dimensional law. Let us first 
recall that, thanks to stationarity and independence of the increments, we have that 


Eexp(-/iSt) = exp [-t4>(p)], t >0, p>0, 


(1.4) 


where 4>(/r) is called the Laplace exponent of the process S^, t > 0. Different expressions for 4>(|U) give rise 
to different subordinators, but not any functional form is allowed. This is exactly what the Levy-Khintchine 
formula for subordinators tells us. 

Theorem 1.1. Any function 4>(/r), /r > 0, that can be put in the unique form 


4>(/r) = bp-i- 


(1 - m(dx), b > 0, 


(1.5) 


where m(dx) is a measure concentrated on (0, oo) such that f (1 A x)m(dx) < oo, is the Laplace exponent of a 

strict subordinator. Conversely, if 4>(/x) is the Laplace exponent of a strict subordinator, there exist a non negative 

number b and a unique measure m with f (1 A x)m(dx) < oo such that (11.51) holds true. 

0 
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For a detailed proof see 0] . The measure m is called the Levy measure and b is the drift associated to the 
strict subordinator S^, t > 0. For the sake of simplicity, in the following we shall use the term subordinator 
meaning in fact a strict subordinator. 

A particularly simple example of a subordinator is the so-called stable subordinator. This is characterized 
by the Laplace exponent a e (0,1), corresponding to a Levy measure m(dx) = [a/r(l — a)] x“^““dx. 

Note that the case a = 1 is omitted as it is trivial. The one-parameter a-stable subordinator is at the basis 
of the probabilistic theory of anomalous diffusion based on subordination. In fact, it is a buildiim block for 
the numerous processes connected to fractional evolution equations and their generalizations 15] and 
also for other processes connected for example to point processes [33, 3]- Moreover, the importance of stable 
subordinators stems also from the fact that they are scaling limits of some totally skewed generalized random 
walks Given a subordinator Sf, t > 0, it is possible to define its right-inverse process [Q] as 


Ef = inf{w > 0: S„ > t}, t > 0. (1.6) 

The process E,- is non-Markovian with non-stationary and dependent increments [|^ |2^ . 

Let us now introduce some basic facts on tempered subordinators, in particular for tempered stable- 
subordinators. First of all, let us refer the reader to the paper for a complete and detailed account 
on the general theory of tempered stable processes. However, it is interesting to note that tempered models 
already appeared in the literature (amongst others, for example, the KoBoL model [I^. Ml). The class of tem¬ 
pered stable subordinators has been introduced in order to have processes possessing nicer properties than 
those of standard stable subordinators. In practice the Levy measure of a stable subordinator is exponentially 
tempered, thus obtaining 


m(dx) = e —rx ^ “dx, a e (0,1), <? > 0. (1.7) 

r(l - a) 

This simple operation produces desirable results as seen immediately by realizing that the Laplace exponent 
in this case can be written as ^(ju) = (i? + For more information regarding tempered stable sub¬ 

ordinators and associated differential equations, see Eb- In Section 121 we introduce a process which has 
both the standard stable subordinator and the tempered stable subordinator as special cases and we study 
its properties. In SectionjS] instead, we describe connections of the introduced process with some difference- 
differential equations involving a generalized fractional difference operator acting in space. An interesting 
special case related to these equations is that regarding the so-called space-fractional Poisson process . 
In the last Section [3(3l we present a further generalization leading to a process which is no longer a Levy 
process. This generalization is based on the study of similar difference-differential equations but involving 
the so-called regularized Prabhakar derivative in time that generalizes the Caputo derivative in time. 


2 A subordinated sum of independent stable subordinators 


Let us consider the filtered probability space (n, P) and the process 




n \ vr ^ 
T] V 


t>0,neN,rj>0,vre(0,l), r = l,...,n, (2.1) 


where t > 0, r = 1,..., n, are n IJ-adapted independent vr-stable subordinators. Let us further consider 
a positive real parameter 5 such that [5] = n and the S^-adapted tempered 5/n-stable subordinator 
t > 0, evaluated at p", i.e. with Laplace exponent 


‘F(m) = (t)"+ - 1?®, 5 > 0, Tj > 0. 

Proposition 2.1. The Laplace transform of t > 0, reads 

Eexp =exp(-t [(t) +/r'')" - tj"] ) , t>0, p>0. 

Proof. Formula (12.3D can be proven by noticing that 

f 


E exp (-p ^= ]~[ E exp 


r=l 


V 


n \ vr 

-Ml .. 1 


( 2 . 2 ) 


(2.3) 


(2.4) 
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and 


by recalling that Eexp /rvf j = exp tp^^. We 


thus obtain 


Eexp (-P.y;-") = flexp r-tp"- 


r=l 




(2.5) 


and hence the claimed result. 


□ 


Let us now consider the probability space and the subordinated filtration 0 = 

5 /n^^^Q. In the following, our intorost will focus on the subordinated und 0 - 3 .d 3 .pt 6 d process 

t > 0, 5 > 0, V e (0, 1), T) > 0. 

Proposition 2.2. The Laplace transform of the process t >0, can be written as 

Eexp =exp(^-t [(rj + p’')® - p®] ] , t>0, p>0. 

Proof. By definition (12.611 we can write that 




( 2 . 6 ) 


(2.7) 


Eexp 


00 

Eexp (-PP e ds] 

Jo 

00 

exp (^-s [(tj + p'')” - T}"]) P e dsj . 

Jo 


( 2 . 8 ) 


Then by using (12.2D we directly arrive at ( 12.7D . □ 

Remark 2.1. Plainly, when 5 = n = 1, the process (12.6D coincides with a standard v-stable subordinator, while 
for V —> 1, 5 & (0, IX it is a tempered 5-stable subordinator. As an example, we also note that for 5 = 2, 
V G ( 0 , 1 / 2 ), = V^'’ + is the rescaled sum of two independent stable subordinators. 

In the following theorem we prove that the ©-adapted process t > 0, is in fact a subordinator and 

derive its associated Levy measure. 

Theorem 2.1. The &-adapted process t >0, is a subordinator. Furthermore, its associated Levy measure 

is 


m(dx) 



(2.9) 


Proof. In order to prove the statement of the theorem we make use of the Levy-Khintchine formula for 
subordinators. The Laplace exponent associated to t > 0, that is '1>(|U) = ( 17 -I-— 17 ^, can be 

retrieved with the following steps. 


( 1 -e-'^^) 


n, r ^ 5 t“(^+n) 

e-" ‘P LrX’" G dx)-;-—dt 

2 n p _ 5^ 


’ r(i+D 

dte ^ - - -— 

r“ r(i+D 

dte-^^ -p-— 

„r(i-f) L 

r°° K i--fi+-) 

' dte-”"'- p -^fl-e-'K”+'''' 

nr(l-f)V J 


(1 _ P ('^■ 7 /v.n G dx) 

D 


dt 


I Jnr(l-f) 


( 2 . 10 ) 
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dt 


0 -' 


n r 

(1-!) 

5 

1 

n r 

5 

(i-H 

1 

n r 

(i-H 


1 f 


'v+i^ 

5 t 

<* 00 

dt(l-e-’’"') 

0 

5 r 


n r ( 

■* 00 

(rj + e-y<^^+>^'y dy r^-^dt - 

-•00 r 


0 J 


0 J 




7]'‘e 


/* 00 /* 00 

[t] dy t~^~ndt — j j r]^e~^'^'"dy t~^~n 

Jo Jy 


dt 


/* 00 




y 

r 00 


n 5 n 

dyr]’'e y " - 
o 


r('-!) 


dz e 


I* 00 


(’i + m'')"; r(i-f) 


dze ^ — 


From the above computation and Theorem ll.il the claimed form for the Levy measure easily follows. □ 

Remark 2.2. Note that the determined measure can also be expressed by means of Riemann-Liouville fractional 
derivatives. Indeed, since formula (2.117)] 


s t-(-U) t-OU) 
„r(i-f)“r(i-(i + f)) 



( 2 . 11 ) 


where 




I* 00 


r(m — j3) dt" 


fiy) 


a (t-j/ 


—m+1 


dy, 13 >0, m=[p] + l, 


is the Riemann-Liouville fractional derivative 13517 . we obtain readily that 


m(dx)= — 


e djc) (^I>o^"lj (t)dt. 


( 2 . 12 ) 


(2.13) 


In order to study the operator associated to the subordinator we start from the Laplace transform 


Ee 




By taking the derivative with respect to time we obtain 

d 
dt 

From Theorem l2.1l we know that 




t >0, p> 0. 


(1 - e-'^^)m(dx), 


[(p + p-f-y] = 
where m(dx) is given by (12.9D . Hence, 

I* 00 

[(17 +/t'')^ - Tj®] Ee“'^’i®'’ = rEe“'^i®'’ - e“'^-*'Ee“'^’>®'' j m(dy), 

Jo 

and equation (I2.15D becomes 


—Ee = — I ^Ee — e '^^Ee j rn(dy). 


(2.14) 


(2.15) 


(2.16) 


(2.17) 


(2.18) 
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Now, a simple inversion of the Laplace transforms leads to 


(* 00 


—G dx)/dx = - G dx)/dx - G d(x - j))/d(x - y)) m(dy). (2.19) 

J 0 

Let us denote now v(x, t) = P g dx^ /dx and with 




V'-'x 


[fix, t)-f[x-y, t)) m(dy) 


the generating form of the operator. The probability density function v(x, t) of satisfies 


—v(x, t) = -^ 0 ^’ v(x, t), x> 0 , t> 0 . 


( 2 . 20 ) 


( 2 . 21 ) 


Remark 2.3. We know [2^ that the Prabakhar derivative is defined, for a suitable set of functions [see\2^ 
for details], as 


(t-[C(t-y)“]/(j)dy, (2.22) 

0 

with 0,17 G C, 01(0) > 0, 01 ( 17 ) > 0, G C, t > 0. The fractional derivative appearing in (I2.22D is the 
Riemann-Liouville fractional derivative with respect to time t and 




'^a,77,C;0+*' 


^ 0 + 


is the generalized Mittag-Leffler function (see e.g. 


^ x''(r)r 

= 1 V K,m,rGC, 01(/c)>O, (2.23) 

’ t!r(jcr + ct) 

The Laplace transform of the Prabakhar derivative is 


f* 00 

""'DLf; 0 +/^Odt = s'’(l-Cs-“)V~(s), s> 0 , 

Jo 

where f is the Laplace transform of f. Therefore, with our choice of parameters, 

J(t)dt = s"S(l + ris-n^f(s) = (s^ + r 7 )V~(s), s > 0 . 
This implies that = D* _^ o+ ~ thus 

I* 00 

e-f^\Ql-^v(x, t)dx = [(m" + r?)" - rj^] v(p, t). 


^v,v5,—ri]0+’^ 


(2.24) 


(2.25) 


(2.26) 


0 


Remark 2.4. Note that an alternative representation for the operator tan be given in terms of Riemann- 

Liouville derivatives. We have 


, 0 r=( 0 +v -0 • 

This can be proved by simply computing the Laplace transform, as follows: 

(* 00 


e-'^" [(p + D-^y-r^^]f(x)dx. 


We can perform a formal expansion by means of Newton’s theorem obtaining 


-M- 2 ( ) p^-’-Dllf(x)dx = [(p.'' + r?)^ - rjS] f(p). 


r=l 


(2.27) 


(2.28) 


(2.29) 


We considered the set of functions for which the semigroup property for the Riemann-Liouville derivative holds, 
and also that 11°^ is the identity function. Note finally that ( I2.29D coincides with (12.261) . 
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3 Generalization of the space-fractional Poisson process 

In this section we p resent a study of possible generalizations of the so-called space-fractional Poisson process 
[fssll (see also [|3^ for other more generalized results and special cases). Furthermore, in the following we 
will see how the process introduced in the previous section arises rather naturally in this framework as a 
time-change of an independent homogeneous Poisson process. Furthermore, Section [331 shows the effect of 
replacing the integer-order time derivative in the governing difference-differential equations with a non-local 
integro-differential operator with a three-parameter Mittag-Leffler function in the kernel, i.e. the so-called 
regularized Prabhakar derivative. 


3.1 Space-fractional Poisson process 

In order to make the paper as self-contained as possible we summarize here the basic construction of the 
space-fractional Poisson process as it was carried out in 113^ . The original aim was to generalize a homoge¬ 
neous Poisson process in a fractional sense by introducing a fractional difference operator in the governing 
equations acting on the state space. The chosen operator is (1 — B)“, a e (0,1) (where B is the backward-shift 
operator, that is, given a function /j. depending on an index fc e Z, B(A) = and S'-(A) = B(B''-i(A))) 
which appears in the study of long memory time series. The introduction of the fractional difference operator 
(1 — B)“ = 0 implies a dependence of the probability of attaining a particular state to those 

of all the states below. The difference-differential equations for the state probabilities of the space-fractional 
Poisson process read 

= Q:e(0,l], A>0, 

[pfc(O) = 

where Pjt(t) = P{JV“(t) = k}, k > 0, t > 0, are the state probabilities of the space-fractional homogeneous 
Poisson process iV“(t), t > 0. If a = 1, the process JV“(t) reduces to the homogeneous Poisson process of rate 
A. The space-fractional Poisson process possesses independent and stationary increments and E(JV“(t))^ = oo, 
h= 1,2,.... The Cauchy problem (13.11) is easily solved using the probability generating function G(u, t) and 
the Laplace transform. It turns out that the probability generating function can be written as 

G(u,t) = e-^“'(i-“5“, |u|<l, (3.2) 


and by means of a simple Taylor expansion the state probability distribution is recognized as a discrete stable 
distribution and reads 


PkM 


fc! 


E 


(-A“t)'' r(ar-fl) 
r! r(ar-|-l —fc)’ 


fc > 0. 


(3.3) 


The behaviour of the process is made apparent by the subordination representation 

Af“(V'/) = JV“’’(t), a, ye (0,1). 


(3.4) 


For a = 1, the time-changed process JV^(y/) increases super-linearly with jumps of any integer size. Following 
the same lines, in the next section we will construct a generalized model by suitably adapting the operator 
acting on the state-space. 


3.2 Generalized model 

Consider the equations 


|^PfcCf) = -{[B + ^''(l--B)'']®-Tj^}p/c(0, vnG(0,l), n= TSl, 5,TjeK+, A > 0, 

[pkW = 


(3.5) 


The generating function G(u, t) = can be easily deduced from (13.51) by considering the following 

steps. 


dt 


f) = -2 ]2 f k-ir 2 u^Pk-mM 


(3.6) 


k=m 


1 




The solution is thus 




r=l 


m=0 


A:=0 


=-E 


vr 


m 


[-urG[u,t) 


-G(u,t)2 Jrj^-'-A-d-ur 

r=l 

-G(u,t) 1 ^( 17 + A''(l-uf)®- 17 ® j , G(u,0) 


= 1. 


G(u 




(3.7) 


The following theorem describes the structure of the point process for which the state probability distribu¬ 
tion satisfies (13.5D . For the proof we will make use of results described in Section[2l 

Theorem 3.1. Let t >0, be the process in (12.61) and N[t\ t > 0, be a homogeneous Poisson process of 

rate A > 0 independent The probabilities ’^) = k}, k > 0, t > 0, of being in state k at time t, 

are solutions to d^ . 

Proof. By exploiting the subordination property, the probability generating function of the process 
can be written as 


e ds) 


Then, by using Proposition [2]^ we obtain the probability generating function 03.7D . 


(3.8) 

□ 


The form of the probability generating function shows that the mean value of the associated process is 
infinite unless the subordinating process reduces to the tempered stable subordinator, i.e. if v = 1, 5 e (0,1). 
From the probability generating function G(u, t), it is possible, by a simple Taylor expansion, to obtain the 


state probability distribution p^it) = P[JV(„5Jf’®) = fc], fc > 0, t > 0. We have the following theorem. 


Theorem 3.2. The state probability distribution Pi^it) = P[JV(„2Jf’^) = fc], fc > 0, t > 0, reads 


00 


r(vm -I-1) 


—I m!r(vm — fc -I-1) 

m =0 ^ ■' 


= pf’i 


fc! 


Z 


r(vm -I-1) 


^ m!r(vm — fc -I-1) 

m =0 ^ ■' 


P ] § r! 
iV't 


(-ti7*)T(r5 -I-1) 


r(r5 — m -I-1) 


(3.9) 


■ (1,5) 


(1 -m,5) 

— rjt 


where ypfjiz) is the generalized Wright function (see 12317 . page 56, formula (1.11.14)). 
Proof. Starting from the probability generating function G(u, t), t > 0, \u\ < 1, we have 

[-t(^r] + r(i-uyfY 


G(u, t) -■ 






Z' 


,07 




00 f .\r 00 

s^z 


m=0 
r 00 




r=0 


m =0 
,h 00 


r=0 
5 r—m T vm 


r! 


(3.10) 


A^Cl - uY 
00 

5 r—mp^vm ^ | 


(-uf 


yA.n'ZPIvZflf’ 

^ h\ ^ rl ^ 


h=0 

r(r5 + l) 


h=0 


Z"“ 

h =0 




h'. ^ r! ^m!r(r5-m+l) 

(- 1 )^^ r(vm-Fl) 

~7ir 


^5r-m^vm r(vm + 1) 


Z 


^ m!r(vm — 7i -I-1) 

m =0 ^ 




r(vm — 7i + 1) 
r!r(r5 —m + 1) f 


and thus the claimed formula 03.9D . 


□ 
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Remark 3.1. For 5 = 1 we obtain the discrete stable distribution (2.15) of characterizing the behaviour of 
the space-fractional Poisson process. Indeed it suffices to compute 


V 


-m trj 


z 

r=0 


(-fqyrfr + 1 ) 

r!r(r — m + 1 ) 




S = (-tr. 


(r — m)\ 


(3.11) 


See also nSdil . Section 4.1. For general information on discrete stable random variables, the reader can refer to 
l32] or 11317 . For v = 1, 5 e (0, 1), a time-change given by a tempered v-stable subordinator is retrieved (see 
[34], Section 4.2). 

Remark 3.2. Since N(t) is a Levy process, the subordinated process is a Levy process with Laplace 

exponent $(/r) = (tj+A'^(1—tj®, and thus it possesses stationary and independent increments. Moreover 
is not in general a renewal process. This comes simply from l247 (see also Hdll ) and from the fact that, 
since is a properly rescaled and time-changed linear combination of independent stable subordinators, its 

inverse process £[ = inf{w > 0 : > ^3} general non-stationary and dependent increments. 


Remark 3.3. The reader can compare the state probabilities (13.9D with the results obtained in 11341 . Section 2, 
with a suitably specialized Bernstein function f. 


3.3 Time-fractional generalization with regularized Prabhakar derivatives 

We saw from the above analysis (see Remark [2)3] and the subordination result) that the difference operator 
(13.5D is in practice connected with the Prabhakar derivative. This suggests that a further generalization, 
which can be still tractable, could involve a regularized Prabhakar derivative acting in time. This is what was 
considered in [17] in the case of the classical difference operator and in [fisll for the generator of a Levy 
process. 

Let us first recall the definition of the regularized Prabhakar derivative. 

Definition 3.1. Let l3,co,Y,a^C, Dl(/3), 31(a) > 0, n = [3l(^)], / e AC''[0, b], 0 < b < oo, where 

( d"-! 1 

AC" [a,b] = if : [a, b] ^ K: ^ n_i / (^) ^ absolutely continuous in [a, b] >. 

The regularized Prabhakar derivative is defined as 


CDt 


+/(x) = D^ 


x'' 




k=0 


(3.12) 


For further details the reader can consult Gi. 

By means of the regularized Prabhakar derivative we can construct the following Cauchy problem. 


CDt 

a.B.—co.O''' 


[Pfc(O) = 5,,o, 


PkW = -{[v + r(^l-ByY-r)^^Pi,(t), vne(0,l), n= [5], 5,r]> 


(3.13) 


with the constraints £u>0, 7 > 0 , 0<a<l, 0</3<l. We also have 0 < ^ [yl/y — ra < 1, V r = 0,..., [y], 
if y 7 ^ 0. With similar calculations as for (13.61) . we obtain for G(u, t) = u^p^.(t) 




a,^,— cl), 0 + 


G(u, t) = -1 [17 -F A’'(l - u)’'] ® - Tj®| G(u, t), vn e (0,1), 


(G(u,0) = l. 

By applying the Laplace transform with respect to time t we have that 

s^-^l-F 


G(u,s) = 


s^(l-Fdus “y-h (r)-h — uyy — r/^ ’ 


which can be written, for [[(17 -F A''(l — u)'')® — ') 7 ^]/[s^(l -F cos “)’"][ < 1, 

CX) 

G(u,s) = 2 {-[(p -F A’'(l - u)’')® - -F ms-“)-"C 


(3.14) 


(3.15) 


(3.16) 


n=0 
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We can now invert term by term the Laplace transform by using result (2.19) of [12211 and Theorem 30.1 of 
Hi . The probability generating function can be written as 


G(u, t) = 2 {- [(r, + - uf- V'] (3-17) 

n=0 

Note that formula (2.20) of and equation ( I3.17D only coincide when y = 0 in ( I3.17D and, at the same 
time, in ( 2 . 20 ) of [fs^ . /(•) is specialized to (17 + A’'(.)'')® — 17 ®. 

By expanding the probability generating function ( I3.17D we can derive the state probability distribution 
for this generalized model. Indeed we have 


IL CAJ 

G(«, t)= i; 


n=0 

Rearranging, we get 


r=0 


m=0 


r5 


m 


h=0 


(-ud. 


r— 00 n 00 . ^ 


n=0 


r=0 


m=0 


r(vm + l) 

^or—mTvm '' 




r(vm — fc + 1) 


(3.18) 


(3.19) 


Remark 3.4. Note that for y = 0, 5 = 1, the generating function ( I3.17D reduces to that of the space-time 
fractional Poisson process 13317 of parameters (/3,v), while for 7 = 0, v = 1, 5 G (0,1), we obtain that of the 
tempered space-time fractional Poisson process. If we set j = 0, 5 = 1, v = 1, we get the time-fractional Poisson 
process of parameter (3; for y = 0, 5 = 1, jS = 1, we obtain the space-fractional Poisson process of parameter v. 

We now show that the probabilities p^ft) of (I3.19D are in fact state probabilities of a suitably time-changed 
homogeneous Poisson process. Consider the stochastic process, given as a sum of subordinated independent 
stable subordinators 


Irl 




r=0 


jjUl- 

r 

$(t) 


t > 0. 


(3.20) 


ji ji 

The random time change is defined as 4>(t) = V/’'', where V/’”' is a stable subordinator independent of 

all the others and where 0</3[y]/y—ra<l holds for each r = 0,1,..., [y]. The hitting time process can 
be defined in turn as 


iff = inf{s > 0 : > t}, t > 0 . 

We are now ready to state the following theorem. 


(3.21) 


Theorem 3.3. Let Ilf, t > 0, be the hitting-time process in formula (I3.21D . Furthermore let be 


a homogeneous Poisson process of parameter A > 0, subordinated by the process defined in 

independent of The time-changed process 


n(t) = iv(^2j;/f), t>0. 


and 


(3.22) 


has state probabilities (I3.19D . 


Proof The claimed result can be proved by writing the probability generating function related to the time- 
changed process n(t) as 


2u¥(n(t) = k)= G dy). 


(3.23) 


fc =0 


Therefore, by taking the Laplace transform with respect to time and taking into consideration Theorem 2.2 of 
we have 


/* 00 


00 


0 J 


e-y[Cri-tAqi-i,r/-r,«]-stp(-^^ e dy)dt = 


sl^-\l + cos-^y 


s^(l -I- cos “)t -I- (17 -I- A''(l — uyy — rj^ ’ 


which coincides with (I3.15D . 


(3.24) 

□ 
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As a byproduct of our analysis we obtain 


00 

=2(-+ t > 0, s > 0, (3.25) 

fc =0 

that generalizes formula (3.37) of [|^ with the Laplace exponent suitably specialized. 
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